Spectral and polarization dependencies of luminescence by hot carriers in graphene 
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The luminescence caused by the interband transitions of hot carriers in graphene is considered 
theoretically. The dependencies of emission in mid- and near-IR spectral regions versus energy and 
concentration of hot carriers are analyzed; they are determined both by an applied electric field and 
a gate voltage. The polarization dependency is determined by the angle between the propagation 
direction and the normal to the graphene sheet. The characteristics of radiation from large-scale- 
area samples of epitaxial graphene and from microstructures of exfoliated graphene are considered. 
The averaged over angles efficiency of emission is also presented. 

PACS numbers: 78.67.Wj, 78.60.Fi, 72.80.Vp 
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I. INTRODUCTION 

Electro- and photo-luminescence of the bulk semicon- 
ductor materials and structures have been used for more 
than fifty years both in the characterization these mate- 
rials and in the operation of light-emitting devices. [l| 
Emission of radiation in the far- and mid-IR spectral re- 
gions caused by nonequilibrium charge carriers also have 
been studied for two-dimensional systems, see reviews in 
Ref. [2( , and for the transitions between the subbands of 
heavy and light holes in p-Ge. [3| In graphene processes 
of emission of radiation are actual due to the effective 
heating of carriers by dc electric field [3-01 and due to 
the the efficient interband transitions excited by photons. 
[7( In this work we analyze spectral and polarization de- 
pendences of emission for interband transitions (on fre- 
quencies exceeding the frequency of relaxation) induced 
by hot carriers in a bipolar graphene. The characteristics 
of radiation are considered for two cases of emission: (a) 
from large-scale-area samples of epitaxial graphene and 
(b) from microstructures of exfoliated graphene, see Figs, 
la and lb, respectively. The analysis is based on the 
quasiclassical kinetic equation for 3D photons where the 
interaction with 2D carriers is described by the bound- 
ary condition at graphene sheet (see Ch. 4 in [8|] and a 
similar approach for the case of acoustic phonon emission 

1). 

Spectral dependencies of emission arc determined by 
the character of distributions of nonequilibrium electrons 
and holes, wherein a heating electric field modifies not 
only an effective temperature of carriers but also electron 
and hole concentrations. These nonequilibrium charac- 
teristics are dependent on a lattice temperature, T, a 
strength of applied electric field, E, and a gate voltage V g . 
Due to direct interband transitions and a linear character 
of the energy spectrum of graphene, a radiation with frc- 
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FIG. 1: Geometry of luminescence by hot carriers from 
graphene sheet of: (a) a large-scale-area sample (S), with 
detector (D) placed in the near-zone; (b) a finite-size sample 
(Si), where G r is the energy flow in the far-zone, (c) Polar- 
ization characteristics of radiation. Here n is the normal to 
2D plane, ei : 2 are the unit vectors of polarization, and Q is 
the wave vector. Angles 6 and cf> define propagation direction 
and in-plane orientation, respectively. 



quency uj is emitted for transitions between electron (e) 
and hole (h) states with the momentum p u = huj/2vw, 
where v\y ~ 10 8 cm/s is the characteristic velocity of the 
Weyl- Wallace model (degenerated both on spin and val- 
ley quantum numbers). [lOj It appears that spectral de- 
pendence of radiation is proportional to the product of 
e- and /i-distribution functions, fep^fhp^, and the max- 
imum energy of emitted photons is determined by both 
the effective temperature of carriers, T c , and the electron 
and hole concentrations. The spectral behavior of radi- 
ation is strongly dependent on the densities of carriers, 
which in turn are essentially dependent on V g . 

The polarization dependence of emission results from 



the chiral character of the neutrino-like states, [Tl| so 
that the matrix elements of interband transition are de- 
pendent on orientation of 2D momentum, p. [5| Due to 
this, a strong dependence of polarization from the an- 
gle 9, between the direction of propagation of radiation 
and the normal vector to the graphene sheet (see Fig. 
1), is obtained. In addition, a weak the current- induced 
anisotropy of distributions, oc {eEr m / 'p) 2 <C 1 (here r m is 
the momentum relaxation time and p is the characteristic 
momentum of carriers) , leads only to a weak anisotropy 
of emission. Such weak variations of polarization are de- 
termined by the angle, <p, with respect to the direction 
of the in-plane applied field, and pertinent contributions 
are omitted below. 

The paper is organized as follows. The basic equa- 
tions governing the emission of radiation in graphene are 
considered in Sec. II. In Sec. Ill, we analyze the polar- 
ization characteristics and the spectral dependencies of 
radiation. The concluding remarks and discussion of the 
assumptions used are given in the last section. 



II. BASIC EQUATIONS 

Emission of radiation due to the interband transitions 
of nonequilibrium carriers in graphene is described by the 
Wigner distribution function of photons, Nq^ with the 
polarization indexes fi = 1,2. Outside of the graphene 
sheet, A*q^ obeys the following kinetic equation [8] 



vq • V T N% = J R 



(1) 



Here vq = duiQ/dQ = cQ/Q is the velocity of photon 
with the frequency u)q, Q = (q, q±) is the wave vector 
of photon with the in-plane (transverse) component, q 
(q±_), and c is the speed of light in the medium around 
the graphene sheet with dielectric permittivity e. The 
collision integral Jr describes the relaxation of photon 
distribution outside of graphene (the glancing photons, 
with q± — > 0, are not considered). The boundary condi- 
tion, around the graphene sheet placed at z = 0, takes 
form 
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where vj_ = cq±/Q and Iq" determines the speed 
of spontaneous emission due to interband transitions. 
Within the collisionless approximation, Iq" is deter- 
mined by the direct interband transitions (/ ^ I') between 
the states \lp) and |/'p), with distribution functions fi p 
and /z'p, and is given by H 
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xflpO-- fi'p)8(fiwQ+ei> p -eip). (3) 

Here the factor 4 comes from the spin and the valley de- 
generacy, L 2 is the normalization area, and the polariza- 
tion vectors, eQ^,, are defined by relations (Q • eQ M ) = 



and (eq^ • e Qfl , 



8^^ . It is convenient to choose 



them with the angles ±7r/4 towards the plane formed 
by the vectors Q and n, see Fig. 1. Point out that in 
Eq. ([3]), due to the energy conservation, contribute only 
the direct interband transitions between the states with 
I = 1 and V = — 1, where the ^-function obtains the form 
S(huj Q - 2v w p)- 

Thus, the characteristics of emission (its intensity, 
spectral dependency, and polarization) are expressed via 

-^Qr determined by Eqs. ([I])-©- The energy flow den- 
sity is defined by the standard formula [3, [l2J 



G r = E/(ga v Q^Q^ 
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The polarization properties are characterized by the 
Stokes parameters, £ x , £ y , and £ z which are introduced 
by relations [12J 
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where Nq = £V Nq 1 determines an intensity of radia- 
tion propagated along Q. Below we restrict ourselves by 
the nonabsorbing medium, Jr — > 0, and consider (a) the 
in-plane homogeneous geometry, which is corresponded 
to a large-area sample of epitaxial graphene, and (b) a 
small-size sample of exfoliated graphene. 

For the case (a), Eqs. (p} and ([2]) have the solution 
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which describes emitted radiation with q± > (q± < 0) 
at z > (z < 0). The only nonvanishing component of 
the energy flow density G± is directed along 0Z and it is 
given by 



G ± = Y. 



J (2^)3 Hu}qI Q = J duJ J dnQ 
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(9_l>0) 

where the differential flow, d 2 G/dujd£l, is introduced 

d 2 G = hiv 3 ^ w 

dujdQ (2tt5) 3 2^ Q ' 



(8) 



In addition, it is convenient to introduce the 

frequency-dependent differential flow dG/duj — 

J dflQ(d 2 G/dujdil) that is averaged over the 

(<?-L>0) 

solid angle J7q subtended by the infinite plane. 

Considering emission from the sample placed within 
the in-plane region S\ [case (b)], we suppose that Nj^ 



is zero outside of Si-region (i.e., we neglect by the edge 
diffraction effects). Thus, one obtains the solution 
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with Mq (x,y) = outside of the Si-region. In the far 
zone, R ^> \fS\ (see Fig. lb), the tangential components 
of G vanish and the radial component of energy flow 
takes form: 
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where the differential flow, introduced by Eq. 
pears. 

It is convenient to rewrite the speed of spontaneous 
emission Eq. ([3]) by making the replacement from the 
band quantum numbers (1,1') to the electron-hole rep- 
resentation, when the distributions fi p are substituted 
for electron and hole distributions as: /i p — > f ep and 
(1 - /-i p ) -> fh-p- Then Eq. ((3J) obtains the form 
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y-fe-pfh-p5(?lUJQ - 2v w p), 
% [tp) = (Ip |e Q/i • v| - lp)* (lp |e Q ^ ■ v| - lp),(ll) 

where the interband matrix elements M^^^tp) are de- 
pendent only on the orientation of the unit vectors of 
polarization and on the angle tp, giving the orientation 
of the momentum, p = (p cos tp,p sin tp).\l^ Neglecting 
a weak anisotropy of distributions f e ,hp, one can use in 
Eq. Ijlip the averaged over the p-plane angle (such an 
average, over the angle tp, is denoted using the overline) 
matrix element 
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where 9 = n, Q, see Fig. 1. As a result, the speed of 
spontaneous emission, Eq. (|11[) , we can rewrite as l£$ = 
/q M , due to its dependence only on u> and 9. Finally, 
completing in Eq. (|11|| the integral over p by using the 
energy cj-function, we obtain the speed of emission as 
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where pe — 2E/n(Uvw) 2 is the density of states and the 
distributions are taken at the characteristic momentum 
p u . Therefore, the spectral and polarization dependen- 
cies of radiation are presented in Eq. (|13p by the separate 
factors. 



III. EMISSION CHARACHERISTICS 

Here we study the polarization or the spectral char- 
acteristics of the luminescence determined by Eqs. 
© or ©, 0, and ([T0]l. respectively. Distribu- 
tions of nonequilibrium electrons and holes are de- 
scribed by the quasiequilibrium Fermi functions fk P ~ 
{exp [(v\yp — jUfe) /T c ] + 1} with effective temperature 
T c and the chemical potentials fik, that determine the 
concentrations of electrons and holes, n e and rih- [5( 
Instead of the concentrations n e ^ it is convenient to in- 
troduce the surface charge eAn — e(n e — nh), that is 
defined by the gate voltage V g , and the total concen- 
tration n — n e + nh, that is defined by a character of 
the generation-recombination processes. In addition, the 
effective temperature T c of the hot electrons and holes 
can be estimated from experimental data |14| and from 
calculations. 0, Q 



A. Polarization of radiation 

First, let us consider the polarization characteristics 
of radiation emitted. Here we use that the spectral and 
polarization dependencies are separated, see Eqs. ([S]) and 
(fiUl) . and the Stokes parameters ([5]) are expressed only 
through M^j/, Eq. (fT"2)) . Therefore the polarization of 
radiation is independent of the frequency or the character 
of the carriers distributions. Then, for the geometry of 
Fig lb, in the far-zone by using Eq. (fT3"| in Eq. ([5]) we 
obtain that £ y = £ 2 = and 
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(14) 



which determines a degree of the linear polarization as 
function of 9. For the geometry of Fig. la for any plane 
wave contribution, with given Q, we again can introduce 
the Stokes parameters Eq. ([5]) and they have the same 
form as for the geometry of Fig lb. 

From Eq. (|14p it follows that for the normal propaga- 
tion, 9 — > 0, the emitted radiation becomes nonpolarized, 
in agreement with the absence of any preferential direc- 
tion over the graphene sheet. If to take into account a 
weak lateral anisotropy of the carriers distributions in- 
duced by the applied electric field E then additional de- 
pendence of the polarization characteristics appears from 
the mutual orientation of the vectors E and Q defined by 
the angle 4>, see Fig. lc. This small addendum also de- 
pends on uj and the form of the distribution functions of 
the carriers. For the glancing propagation, at 9 — > 7r/2, 
Eq. (fill shows that emitted radiation is fully linearly po- 
larized, in a 2D plane parallel to the graphene sheet. Ap- 
pearance of the universal angular dependence Eq. (|14p 
allows for separation of the interband contributions in 
graphene from any possible background emission (e.g., 
from a substrate, a cover layer or a gate). 



B. Spectral dependencies 

Next, we study the spectral dependencies and the effi- 
ciency of radiation emitted by hot carriers for the geome- 
tries shown in Figs, la and lb. Performing the summa- 
tion over polarization of Eq. (|T3|) . we present the differ- 
ential flows dG/dtu and d 2 G/dujd£l, see Eq. ©, through 
^ 7^. These differential flows give the relevant inten- 
sities of radiation, by Eqs. and flT0|) . The differential 
flow Eq. ([5]) obtains the form 
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where it is introduced the characteristic density of en- 
ergy G c = y/ee 2 (T c /hc) 3 /n 2 and the spectral behavior is 
determined by the following dimensionless function 
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The angular dependence in Eq. (I15[) is given by the factor 
(1 + cos 2 9); i. e., the intensity of the glancing emission 
there is two times smaller than the normal one. The 
differential flow Eq. (|T5|) integrated over the solid angle 
of half-space gives 



dG 8vr „ 
— = —G C A 
duj 3 



TlLU 
2T r . 



(17) 



i.e., dGjduj is expressed through the function Eq. (ITBl . 
For T c =300 K and e ~3 we estimate the characteristic 
density of energy, G c oc T 3 , as G c ~ 1.06 x 10~ 17 J/cm 2 . 



maximum of radiation, localized at Hlo ~ T c , is slowly 
shifted. For high frequencies A(Tiw/2T c ), Eq. (jTHJ), is 
exponentially decreasing as exp(— Huj/T c ) and at low fre- 
quencies this function grows oc w 3 (the latter asymptotic 
dependence is shown in Fig. 2 by the dashed curve, for 
n/riT a = 1.5). 
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FIG. 3: (Color online) Dimensionless function A versus 
Twj/2T c for doped graphene at T c =300 K and at differ- 
ent generation-recombination conditions determined by ratio: 
n/An =1.25 (1), 2 (2) and 3 (3). Doping levels are governed 
by gate voltages: V a =5 V (a), 10 V (b), and 20 V (c). 




FIG. 2: (Color online) Dimensionless function A versus 
hu/2T c for the intrinsic graphene case at different generation- 
recombination levels determined by ratio: n/nr c =0.5 (1), 1 
(2) and 1.5 (3). 

Thus, the spectral characteristics of radiation are ex- 
pressed through the function A(Tiui /2T C ) and it is shown 
in Fig. 2 for an intrinsic graphene. Here, one defines the 
degree of noncquilibrium of the electron-hole pairs con- 
centration as the ratio, n/riT c , of the total concentration 
n to pertinent equilibrium one (at T c and /j. e = fih = 0), 
where tit c — (Tr/3)(T c /hvw) 2 - This ratio also charac- 
terizes the effectiveness of generation-recombination pro- 
cesses. It is seen from Fig. 2 that as n/riT grows the in- 
tensity of emission essentially increases and the spectral 



For the doped graphene the spectral maximum is 
shifted to higher energies due to the Pauli blocking effect 
arising as the gate voltage, V g , increases (compare Figs. 
3a, 3b, and 3c plotted for T c =300 K). Here the total 
concentration of carriers is defined by the ratio n/An, 
where eAn gives the surface charge density controlled 
by V g . Present calculations are conducted for a typical 
graphene structure on SiC>2 substrate of the thickness 300 
nm. As in intrinsic graphene, for growing concentration 
the intensity of emission increases, moreover, for growing 
V g the maximum also rapidly becomes larger. For high 
energies the spectrum of emission becomes exponentially 
decreasing. 

For increasing temperatures the character of depen- 
dences from n/An and V g is not modified, see Figs. 4a-c 
plotted for T c =600 K. Here, for two times larger tem- 
perature the altitude of A becomes about one order of 
magnitude smaller. However, as G c oc T 3 , the maximal 
differential flows Eqs. (fT5j) and (|17p arc wcakcly modified 
for such changes of T c . In addition, point out that here 
the position of maximum is more strongly dependent on 
the character of generation-recombination processes than 
for the intrinsic graphene (compare Fig. 2 with Figs. 3 
and 4). 
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FIG. 4: (Color online) The same as in Fig. 3 at T c =600 K 
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FIG. 5: (Color online) Energy flow density G± versus 
V g for temperature T c =300 K at at different generation- 
recombination levels determined by n/An =1.25 (1), 2 (2) 
and 3 (3). The dashed curves show the parabolic approxima- 
tion, G±_ OC Vg. 

flow in the far-zone Eq. (fl 
35 1 ! 



G r {6) 
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is expressed through G± as 
G± (1 + cos 2 (9) . (19) 



C. Energy flow 



Here the angular dependence coincides with that of Eq. 
(IT51) and the energy flow (|T9"f decreases oc R~ 2 . Due to 
this, on macroscopic distances, for ^JSi/R 2 ~ 10~ 3 , it 
is possible to register the energy flow ~ 1/iW/cm 2 , while 
G± ~ lW/cm 2 is rather easily achievable according to 
the above treatment. 



Now, let us consider the integral energy flow for the 
cases (a) and (b), when Eqs. ([7]) and ([TO)) can be ex- 
pressed through the differential flow, Eq. (fL3|) . For the 
homogeneous case (a) the energy flow density takes form: 
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For intrinsic graphene at T c =300 K we calculate from 
Eq. (HHJ) that the flow G± ~0.48 mW/cm 2 , 2.02 
mW/cm 2 , and 4.69 mW/cm 2 corresponds, respectively, 
to n/riT =0.5, 1, and 1.5; these n/nx are used in Fig. 2. 
For T c =600 K and the same values of n/nr, we obtain 
that the flow Gj_ ~ 7.68 mW/cm 2 , 32.3 mW/cm 2 , and 
75.0 mW/cm 2 . 

Fig. 5 shows that for doped graphene G± increases for 
growing concentrations of carriers, defined by the gate 
voltage, whereas n/An (or generation-recombination 
level) is fixed. From Fig. 5 it is seen that for V g ^20 
V, which corresponds to the concentration of carriers 
> 2 x 10 12 cm 2 , the energy flow reaches the values 
~ lW/cm 2 . In addition, the dependences G± versus 
V g are close to parabolic ones (see the dashed curves in 
Fig. 5). The temperature dependence of G± is weak 
(less than 10% for the T c increasing from 300 K to 600 
K; pertinent graphs are not shown) because the spectral 
dependences are determined by dimensionless parameter 
huj/T c . 

For the geometry (b) the radial component of energy 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

Summazing the consideration performed, the examina- 
tion of luminescence caused by the interband transitions 
of hot electrons is presented. It is found, that the uni- 
versal frequency-independent polarization of emission is 
realized for a weakly anisotropic distributions of nonequi- 
librium carriers. The spectral dependences of radiation 
are determined by the factor oc uj 3 multiplied by the 
product of electron and hole distributions. These data, 
together with measurements of the angular dependences 
and the integral intensity of radiation (the latter strongly 
depends on generation-recombination processes), allows 
the effective characterization of hot carriers. In addi- 
tion, obtained efficient emission of hot carriers in mid-IR 
spectral regions opens a possibility for using the electro- 
luminescence of graphene as a source of radiation. 

Now let us discuss recent experiments on emission by 
hot carriers from the back-gated transistor structures, 
subjected to a strong in-plane electric field, [14j where 
measured spectral dependences are interpreted by us- 
ing the Planck's law. For the latter the high frequency 
asymptotic coincides with the asymptotic of Eq. (|16p . 
for the case of the quasiequilibrium Fermi distributions. 
From these spectral dependences in [14j it was found the 
relation of T c with a power of the Joule heating. However, 
the character of recombination, that could be determined 
from the intensity of radiation, was not investigated; de- 
pendences on a size of sample, angular dependences and 



polarization characteristics of the radiation also need a 
special treatment. 

Next, we list and discuss the assumptions used. In 
the study of interaction of a radiation with the graphene 
sheet the only simplification used is the neglect by atten- 
uation of radiation propagated along the layer, so that 
obtained results are not applicable for — > tt/2. More 
rigid constraint is imposed by modeling of the distri- 
butions of nonequilibrium carriers via the quasiequilib- 
rium Fermi functions, with given effective temperature 
and concentrations of carriers. This approximation shows 
that important information on the mechanisms of energy 
relaxation and recombination can be obtained from the 
spectra of luminescence. However, for more precise cal- 
culation of the spectral dependences more realistic distri- 
bution functions of carriers will be needed. Further, the 
spectral and the polarization dependences can be sepa- 
rated only in an approximation of a weak anisotropy of 
the distribution of carriers. But such a separation can 



be broken under a strong enough electric field when an 
essential anisotropy of the distribution functions appears. 

We also have restricted our consideration by the ho- 
mogeneous geometry (without taking into account of 
the edge effects) and the far-zone region geometry [ap- 
proaches (a) and (b) in Figs, fa and lb], however, the 
general treatment implies solving of Eqs([T]), © f° r a spe- 
cific geometry that is not well enough approximated by 
any of these two limit geometries. 

To conclude, obtained results show that spectral, angu- 
lar, and polarization dependences of the electrolumines- 
cence provide a convenient method of characterization of 
the hot carriers in graphene (along with the electrooptical 
measurements [15| and a study of the Raman scattering 
[lj, Ha]). Therefore present results will stimulate sub- 
siquent experiments and their theoretical interpretations 
designated for a verification of the relaxation mechanisms 
of nonequilibrium carriers under their heating both for 
the electric field and for the interband photoexcitation. 
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